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Adler function from i? e e (s) measurements: experiments vs QCD theory. 

A.L. Kataev a * t 

a Institutc for Nuclear Research of the Academy of Sciences of Russia, 117312 Moscow, Russia 

An experimentally motivated QCD analysis of the behaviour of the Adler D-function in the Euclidian region 
is described. It is stressed that by taking account of fe-quark mass-dependent o^-effects one obtains better 
agreement between theoretical predictions and experimentally motivated behaviour of the D-function at large 
Euclidian momentum transfer. A more detailed analysis of QCD predictions, including information on quark and 
gluon condensates, requires more precise data on e + e~ — > hadrons, particularly in the energy regions E < Mj/ 9 
and Mj/y < E < 3.6 GeV. Use of experimental determination of the D-function to test the generalized Crewther 
relation is outlined. 



1. Introduction. The process e + e~- 

annihilation into hadrons is one of the most infor- 
mative processes in elementary particle physics. 
Over the past few decades special attention has 
been paid to detailed theoretical and experimen- 
tal study of its basic characteristics, and in par- 



ticular, of the ratio R e 



hadrons) /a(e + e 



H+u- 



(s) = a to t(e + e — ► 
). As a result, impor- 



tant information about the properties of hadrons 
and of their constituents (quarks and gluons) was 
obtained. For example, e + e~ collisions enabled 
discovery in Novosibirsk of the first direct evi- 
dence for the p°-meson [ m. The J/\l/ resonance 
[ 0] was simultaneously discovered by the e + e~- 
collider SPEAR at SLAC, and at BNL as a result 
of study of the process p + p — * e + e~ + X [ ||. 
Rapid confirmation came from a slight increase 
of the beam energy of the ADONE e + e~ collider 

II- 
The observation of \E r and \& particles and 

their interpretation as cc bound states of quarks 
(e.g. see the review [ ||) has served as an incite- 
ment to the development of the theory of strong 
interactions and QCD in particular. It brought 
up the question of the possibility of finding the 
signal for the fifth quark in e + e~-annihilation. 
Moreover, an indirect determination of the mass 
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of the sixth t-quark was made from detailed fits of 
the experimental data of the high-energy e + e~- 
collidcrs LEP and SLC at the Z°-pole taking into 
account the effects of its virtual propagation (for 
a review see [ g]). Impressively, the extracted 
mass of the t-quark turned out to be in agree- 
ment with the its direct "measurement", result- 
ing from subsequent discovery of the 6th quark at 
the Tevatron. 

Among present theoretical studies of data from 
e + e~-colliders data are attempts to find the phe- 
nomcnologically allowed window for the mass of 
the Standard Model Higgs boson. These stud- 
ies are based on combined fits of data from the 
high-energy colliders LEP and SLC (for a re- 
view see [ M) taking into account the effect of 
running of the inverse electromagnetic coupling 
constant from its low-energy value a _1 (0) ~ 
137.0356... to Mz- At this reference scale its 
value was determined using a compilation of the 
available e + e~-data [ g). The result of Ref. [ [|] 
oT x (M z ) = 128.896(90) was recently updated to 
give ar\Mz) = 128.913(35) [§ (for a discussion 
of the fast developing situation in this area see the 
review of Ref. [ [To| and other related works on the 
subject [ 0). Further estimates of experimen- 
tal and theoretical uncertainties in this important 
quantity are on the agenda. 

2. Electron-positron annihilation: exper- 
iment vs QCD predictions. From the per- 
spective of the scanned energy regions, and the 
physical problems under investigation, all past, 



Table 1 

Classification of constructed and planned e + e _ -colliding machines 



Group Machine 



Location Start Data 



Status in 1999 
June 



Beam Energy (b.e.) Region 
E b . e . + E b . e . [GeV] 



AdA 


Frascati 


1960 








ORSAY 


1961 


- 


0.2 + 0.2 


VEPP-2 


Novosibirsk 


1966 


- 


0.2-0.55 + 0.2-0.55 


I ACO 


ORSAY 


1966 


- 


0.2-0.55 + 0.2-0.55 


ADONE 


Frascati 


1969 


- 


0.7-1.55 + 0.7-1.55 


VEPP-2M 


Novosibirsk 


1974 


working 


0.2-0.67 + 0.2-0.67 


DCI 


ORSAY 


1976 


- 


0.5-1.7 +0.5-1.77 


(^-factory 










DA$NE 


Frascati 


1999 


launched 


0.51+ 0.51 


CEA 


Cambridge (USA) 


1971 


- 


1.5-3.5 + 1.5-3.5 


SPEAR 


SLAC 


1972 


- 


1.2-4.2 + 1.2-4.2 


DORIS 


Hamburg 


1974 


- 


1-5.1 + 1-5.1 


II VEPP-4 


Novosibirsk 


1979 


planned to restart 


1.5-5 + 1.5-5 


CESR 


Cornell 


1979 


working 


3-8 + 3-8 


BEPC 


Bejing 


1991 


working 


1-2.5 + 1-2.5 


c-t- factory 


Bejing 


>1999 


was planned 


1.5-3 + 1.5-3 





B-factory 






asymmetric colliders 




III 


PEP-II 


SLAC 


1999 


launched 


9 + 3.1 




B-factory 


KEK 


1999 


launched 


8 + 3.5 




PETRA 


Hamburg 


1978 


- 


5-19 + 5-19 


IV 


PEP 


Stanford 


1980 


- 


5-18 + 5-18 




TRISTAN 


KEK 


1987 


- 


25-30 + 25-30 




SLC 


SLAC 


1989 


working 


45-50 + 45-50 


V 


LEP-I 


CERN 


1989 


- 


45-47 + 45-47 




LEP-II 


CERN 


1995 


working 


65-100 + 65 -100 




TESLA 


DESY 


> 2005 


linear colliders 


250-500 + 250-500 


VI 


NLC 


SLAC 


> 2005 


( projects under 


250-500 + 250-500 




JLC 


KEK 


> 2005 


discussions) 


250-500 + 250-500 



present and proposed e + e~-colliders may be di- 
vided into six groups. We present this classifica- 
tion in Table 1, updating the material given in 
Ref. [ [12]. The new information was gained, in 
part, from material in Ref. [ |L3| . 

These colliders provide complementary infor- 
mation about the behaviour of the e + e~ — ► 
hadrons total cross-section, at different energies, 
from machines of different luminosity. More- 
over, some important physical problems under 
investigation necessitate more precise experimen- 
tal data, not only at high energies, but also in 
low and intermediate energy regions. Indeed, the 
study of the latter regions, both by experimen- 



tal and theoretical methods, may provide better 
estimates of 



acq 



Aa had (q 2 ) = -^-Re I ds- 

•JK Jim, 2 s \ s — Q le ) 



R e+e~ 



(s) 



(1) 



which is the hadronic vacuum-polarization contri- 
bution to the value of the effective fine structure 
constant. 

Another important characteristic of the process 
e + e _ -hadrons was introduced in Ref. [ M . It con- 
cerns the Euclidian Adler D-function 



D(Q 2 ) = Q 2 



R e 



(«) 



w (s + Q 2 ) 2 



(2) 



which can be related to Eq.(l) as follows: 

. . 3tt d , . 

D(-s) = — s— Aa had {s) 
a as 



(3) 



The study of this quantity has a number of at- 
tractive features. Its behaviour was first analyzed 
from an experimental point of view some time ago 
[ II, H) . In the recent work of Ref . [ |§ this prob- 
lem was reconsidered using a compilation of e + e~ 
experimental data [ ||j. The results are depicted 
in Fig.l and Fig. 2 where the shaded areas repre- 
sent the ±lcr band obtained from the data. 
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Figure 1. The "experimental" curve for the Adler 
function together with QCD predictions 
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Figure 2. The low energy "experimental" and 
QCD behaviour of the Adler function 
It is interesting to compare the "experimental" 
behaviour for the D-function with its QCD theo- 
retical expression. One can express it in the fol- 
lowing form 



D^ CD {Q 2 ) = D PT (Q 2 ) + D NP (Q 2 ) 



(4) 



where the first part describes the perturbative 
QCD contributions, while the second term takes 



3 It would be interesting to include in this and related anal- 
yses the new measurements of R e e (s) at y^i values of 
2.6, 3.2, 3.4, 3.55, 4.6 and 5 GeV, from BEPC [Q. 



account of higher- twist effects, which are related 
to the vacuum condensates of quark and gluon 
fields [ [ll| . In the work of Ref. [ [l6| we considered 
an a 2 , mass-dependent expression for D PT (Q 2 ), 
namely 

D PT (Q 2 ) = D^(Q 2 )+D^(Q 2 )^^1 (5) 



D (2) {Q 2 ){ ^m )2 



The Born term expression is well known and has 
the following form 



D^(Q 2 ) = 



3 H < 3/( 1_{ " 

/ 

12x 2 



(6) 



Vl + 4x 2 



m 



VT + 4E + 1 



s/T+4x-l 

where x = m 2 f /Q 2 , while rrif is the pole quark 
mass. The analytic mass-dependent expression 
for DW can be obtained from the calculations 
of Ref. [ nffl , which are in agreement with the re- 
sults of Rcf.[ pOl. Since the mass dependence of 
the 3-loop term D^-term in Eq.(5) is not yet 
known analytically, detailed information about its 
small-mass [ EU E3] and heavy-mass [ E3| expan- 
sions was used in Ref. [ uM to approximate the be- 
haviour of D' 2 ' in the Euclidian region, with high 
precision, using conformal mapping and Pade im- 
provement [ |24|. However, the theoretical expres- 
sion for mass dependence of the double-bubble 
three-loop vacuum polarization diagram was not 
included in the theoretical expression for D^>- 
term in Ref. [ [lq] . There are arguments, based on 
the calculations of (m 2 /<5 2 )a 2 -corrections [ gjj , 
that this contribution is small. Let us consider 
the non-perturbative contributions to D NP (Q 2 ). 
In Ref. [ |l6| the following expression was used 

D NP {Q 2 ) = 3]TQ 2 (8^) 
/ 

1 / n a s ( p 2) \ < ^G 2 > 
12^ 18 i j Q 2 

Otsip 2 ) 



(7) 



2 1 



3tt 

4 a s {^ 2 ) 

27 7T 



+ 0(ai) 



2 ^ ^ < m q qq > 



Q 4 



^))E^^ 



Table 2 
The origin of uncertainties for the "experimental" Adler function at Q 



2.5 GeV and M z 





L>(2.5 GeV) 


rel. err. 


abs. err. 


D(M Z ) 


rel. err. 


abs. err 


Resonances 


.688 (.025) 


3.6% 


0.8% 


.004 (.000) 


5.2% 


0.0% 


E < Mj/n, 


1.068 (.127) 


11.9 % 


4.2% 


.002 (.000) 


14.9 % 


0.0% 


Mj/y < E < 3.6 GeV 


.178 (.035) 


19.9 % 


1.2% 


.001 (.000) 


19.8 % 


0.0% 


3.6 GeV < E < M T 


.850 (.055) 


6.4% 


1.8% 


.032 (.002) 


7.0% 


0.1 % 


My <E< 12 GeV 


.088 (.008) 


8.7% 


0.3% 


.024 (.002) 


9.0% 


0.1 % 


E < 12 GeV data 


2.871 (.146) 


5.1% 


4.8% 


.063 (.003) 


4.9% 


0.1 % 


12 GeV < E QCD 


.162 (.001) 


0.3% 


0.0% 


3.755 (.006) 


0.2% 


0.2% 


total 


3.033 (.146) 


4.8% 


4.8% 


3.818 (.007) 


0.2% 


0.2% 



The terms beyond leading order in a s , calcu- 
lated in the works of Rcf. [ M, turned out to 
be not so important at the moment, because 
the current experimental data for the £>-function 
gives only rough constraints on the contributions 
of dimension-four condensates, which were var- 
ied in Ref.[ 161 within the following intervals 
< ^-G 2 >w (0.336-0.442 GeV) 4 , < m u uu > = < 
m d dd >= -(0.086 - 0.111 GeV) 4 , < m s ss >= 
-(0.192-0.245 GeV) 4 . Nevertheless, in Fig.l and 
Fig. 2 one can see the region, where the addition 
of the nonperturbative QCD corrections to the 
three-loop perturbative QCD expression leads to 
deviation from the experimentally allowed region 
for the D-function at low Euclidian momentum 
transfer, Q 2 . It is possible, that taking into ac- 
count higher-order power corrections [ Eq| might 
improve the agreement with "experimental data" 
at low energies, shown in Fig. 2. Moreover, it is of 
real interest to update the previous analysis, in 
Rcfs.[ £7| Pq| , of the low-energy e + e~ experimen- 
tal data ( see Refs.[ £7| E8|) with the help of the 
QCD Borel sum rule method [ [if] 



R th (s)e- S / M2 ds 



R' 



~{s)e- s ' M2 d S 



In the process of such an analysis, new low-energy 
experimental data, obtained in part in Novosi- 
birsk, can be used. In the theoretical part of 
sum rule of Eq.(8) one should also include per- 
turbative contributions to the coefficients func- 
tion of quark and gluon condensates of dimen- 
sion four [ |2^] , available from the results of Ref . [ 
Bq |, higher-dimension condensates, and informa- 
tion about massless a% contributions to R th (s) 



[ E9|. As to the sum-rule analysis of the high- 
energy e + e~ experimental data, one can try to 
update the studies of Ref. [ BOl, performed with 
the help of the finite energy sum rules approach 

/ R th {s)ds= / R e e ~(s)ds . (9) 

Jo Ja 

Concerning the curves for the D-function ex- 
tracted from experimental data (see Figs. 1,2 of 
Ref. [ [lq]), it is worth noting that the theoreti- 
cal analyses of the recent works [ pU 32 result 
in a description of the low-energy tail of Fig. 2 
by completely different ideas, related to the con- 
cept of "freezing" of the QCD coupling constant 
at small energies [ 33, |34|]. These ideas have re- 
alizations, distinct from conclusions in Rcf. [ |3q| , 
based on application of the PMS approach [ |36fl 
to the four-loop massless theoretical expression 
for R e e (s). Indeed, it was argued in Rcf.[ |3^ |, 
that the observed "perturbative" freezing may be 
spurious, indicating breakdown of a next-to-next- 
to-leading PMS expression, around the scale of p- 
meson mass. Note, however, that the low-energy 
results for the fit of e + e~ data on R e e (s), per- 
formed in Ref. [ [35[ with the help of the following 
■( 8 ) sum rule [M 



R(s,A) = - 

IT 



R< 



(*) 



rds 



(10) 



/o (s'-s) 2 + A 2 
merited serious attention. It should be stressed 
that the low-energy Novosibirsk data of Ref. [ |3!| 
turned out to be essential in this analysis. For 
example, they served as an ingredient in the 
phcnomenological part of the considerations of 
Rcf. [ BO], devoted to the analysis of "analyti- 
cal" freezing of the QCD coupling constant a s 



in the Minkowskian region, and also in Ref. [ [4l| , 
devoted to consideration of the Crewther rela- 
tion [ J42| and its A/5-scheme QCD generaliza- 
tion [|43[, using commensurate scale relations [ H 
(for reviews see [ [45l 46 ). It should be recalled 



that the a^-generalization of the Crewther rela- 
tion connects a massless a^. theoretical expression 
for the e + e _ -annihilation D-function [ E9[ with 
the massless theoretical expressions for the Gross- 
Llewellyn Smith sum rule of the vN deep-inelastic 
scattering and the Bjorken sum rule of polarized 
deep-inelastic scattering, which were calculated 
at the a^-level in Ref. [ E7|. This connection in- 
volves the first and second terms of the two-loop 
approximation to the QCD /3-function. We will 
consider this problem in more detail in the next 
Section. 

We now return to the results of Ref. [ [L6| . Be- 
fore discussing the comparison of the three-loop 
massive theoretical QCD predictions for the D- 
function with the "experimental" behaviour of 
the Z?-function at higher momentum transfer, it 
is worth mentioning that in the process of ex- 
tracting the "experimental" behaviour of the D- 
function not all e + e~ data up to 40 GeV were 
applied. Indeed, in Ref.[ [16) a conservative atti- 
tude was adopted: the real data were replaced by 
perturbative QCD results in certain regions, but 
only where it was obviously safe to do so. Thus, 
in the regions from 4.5 GeV to My and above 
12 GeV perturbative QCD results were used, in- 
cluding massive three-loop [ E2] and a massless 
four-loop QCD contribution [129]. 

The origins of uncertainty in the "experimen- 
tal" curves for the D-function were analyzed 
in Ref. [ Q (see Table 2). The main sources 



come from the region E < M 



•//*! 



accessible 



for more detailed experimental inspection 
at VEPP-2M and BEPC, and the region 
Mj/y < E < 3.6 GeV, which is the priv- 
ilege of BEPC, VEPP-4 and a possible fu- 
ture c — t factory However, as was shown in 
Ref.r 



16 



even at the current level of experimen- 
tal precision one can already obtain new infor- 
mation, namely a demonstration of the impor- 
tance of two-loop heavy-quark mass-dependent 
effects for the "experimental" behaviour of the 
D-function at moderate and high Euclidian mo- 



mentum transfers. Indeed, after including mass 
effects, both in the three-loop perturbative part 
of the D-function, and also in the two-loop run- 
ning of the QCD coupling constant a s from the 
value of a s jjg{M z ) = 0.120 ± 0.003 to lower en- 
ergy scales, via a variant in Ref.[ 48 of the mo- 
mentum n subtraction scheme, one can observe 
the appearance of real agreement of the three- 
loop massive theoretical expression for the D- 
function with the experimentally-motivated Eu- 
clidian curves of Fig.l and Fig. 2. It should be 
stressed, that if we had not included the three- 
loop massive term, a discrepancy with the "exper- 
imental" behaviour of the D-function might have 
been interpreted as requiring non-perturbative 
power corrections from Eq.(7). While the addi- 
tion of the twist-4 power corrections can be of 
real importance in the region of small enough Q 2 , 
deviation of the two-loop curves from the "ex- 
perimental" corridor for the D-function at high 
momentum transfer (see Fig.l) can be associated 
with omission of perturbative QCD contributions. 
Another interesting observation is that the curves 
of Fig.l and Fig. 2 turned out to be rather smooth 



16 1, lacking the resonance enhancements and 



"threshold steps" , typical of the Minkowskian re- 
gion. In view of this we think that possible future 
applications and improvements of the results ob- 
tained in Ref. [ [16| can be useful for more detailed 
tests of perturbative and non-perturbative theo- 
retical QCD predictions in the Euclidian region. 
One such application is presented in the next Sec- 
tion. 

3. New tests of the generalized Crewther 
relation. The Crewther relation [ 42 con- 
nects the amplitude of ir° — > 77 decay with 
the product of the e + e~ annihilation D-function 
and deep- inelastic scattering sum rules, namely 
with Bjorken sum rule of polarized lepton-nucleon 
scattering 



Bjp(Q 2 ) 



6 gv 



Cbjp(Q 



(11) 



[g?(x,Q 2 )-gt n (x,Q 2 )]da 







4 For earlier discussion of the advantages of MOM-typc 
schemes, taking account of threshold effects, see Ref. [W9| 



or with the Gross-Llewellyn Smith sum rule of 
vN deep-inelastic scattering 



GLS{Q 2 ) = 3C GLS (Q 2 ) = / F 3 (x,Q< 

Jo 



)dx (12) 



where the analytical massless perturbative the- 
ory expression for the D(Q 2 ) -function is known 
at the af-level [ £9| , while the massless analytical 
perturbative theory expressions for Cbj p (Q 2 ) and 
Cgls(Q 2 ) up to a^-corrections are known from 
the calculations of Ref. [ [47J . It should be stressed 
that all the quantities we will be interested in are 
defined in the Euclidian region and all the results 
were obtained in the M S'-scheme. In this scheme 
the o^ generalization of the Crewther relation, 
discovered in Ref.[ O, has the following form: 

C B]P (MQ 2 ))C D (a s (Q 2 )) = l+ ^ (Q,) P(a,)(13) 

where Cd{® s {Q 2 )) is coefficient function for the 
Adler D-function, normalized to unity at lowest 
order, j3^ 2 '(a s ) is the two-loop approximation of 
the QCD /^-function and P(a s ) is a polynomial 
starts from a s and contains two terms. It was 
suggested in Refs.[^, [50] that the factorization of 
the /3-function will persist generally, to all pertur- 
bative orders, and can be related to the effects of 
violation of conformal symmetry by the renormal- 
ization of massless QCD. This was later proved in 
Ref. [ pTl . The theoretical properties of the gener- 
alized relation, written down in the M S'-scheme, 
were discussed in detail in Refs.[|4|, [lj| [H] and 
we will avoid their description in this work. How- 
ever, we will concentrate on some phenomenolog- 
ical applications. It is rather useful to use for this 
purpose commensurate scale relations [ [44J , which 
combine the concept of effective charge s [ p2|| (or 
scheme-invariant perturbation theory 
variants of the BLM approach 



153] ) with 
|55], allowing 
one to write the following generalization of the 
Crewther relation [ |d| in the region where the 
heavy-quark masses can be neglected: 



1 



i 

3£^f 



^D(Q 2 )C GLS {Ql) « 1 



(14) 



D(Q 2 )C Bjp (Ql) « 1 



where 










, Ql 


= ~l + ^ 


(15) 




! fa GLS (Ql)\ 


'11 56 
( 12 + T C3 ~ 


- i6Cf)A) 




56 808 
"27"~ C3 + 


320/ 







As specimen physical input we will first use 
the recently obtained value of the Gross-Llewellyn 
Smith sum rule [ p6| 



GLS(12.59 GeV 2 ) = 2.80 ± 0.13 ± 0.17 



(16) 



which does not contradict the less precise similar 
results, obtained using the extrapolation of the 
data [ [57) . Using the results of Ref. [ 58 we con- 
clude that, in this energy region, / = 4 numbers 
of flavours are manifesting themselves. Moreover, 
at momentum transfer Q 2 = 12.59 GeV 2 we will 
neglect c-quark mass effects, which are suppressed 
by a factor m 2 /Ql < 0.19. Thus we conclude that 
the value of the corresponding effective charge is 



&gls{QI 



1 



GLS{Ql) 



(17) 



= 0.067 ±0.043 ±0.06 



where the first (second) error is related to the sta- 
tistical (systematical) uncertainty of the experi- 
mental result of Eq.(16). Using Eqs.(14)-(17) one 
can get the following estimate 



D(Q« 1.8 GeV) w 3.57 ±0.3 



(18) 



which crosses the upper part of "experimentally" 
motivated curve of Fig. 2 at this reference scale. 
A similar conclusion emerges in the case of anal- 
ogous treatment of the experimental result of the 
SMC collaboration for the polarized Bjorken sum 
rule [ p9(l , which is 



Bjp(Q 2 = 10 GeV 2 ) 



0.195 ±0.029 
6 l 9v 



(19) 



1 - 



aBjpjQl) 



As a result of application of discussions in Ref.[ 
|6Q | and considerations presented above, we obtain 
the following estimate for the D-function 



D(Q w 1.6 GeV) « 3.59 



-0.45 



(20) 



which is in rough agreement with the "experimen- 
tally" motivated Euclidian results of Ref.[ fig , 
though with larger errors. 

Of course, our considerations are not so de- 
tailed as the ones given in Ref.[ [51], where an 
expression for the value of the coupling constant 
a s was extracted from the Bjorken sum rule data 
(for earlier discussions of this problem see Ref. [ 
p2[). However, we believe that even approximate 
tests of the generalized Crewther relation, written 
down in the form of commensurate scale relations 
of Ref. [ f41| directly in the Euclidian region, 
give one the feeling that estimates obtained from 
deep-inelastic sum rules are less precise than di- 
rect extraction of the behaviour of the D-function 
in the Euclidian region. This might suggest a 
need for adding higher-twist contributions to the 
perturbative generalization of the Crewther rela- 
tion. 
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